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Abstract 

— H ■ 
> . 

■ It is shown that the rare decays Z uuj and Z — > uuj^ are useful to put model- 

independent bounds on neutrino-one-photon and neutrino-two-photon interactions. 
The results are then used to constrain the r neutrino magnetic moment fi^^ and the 
' double radiative decay uj — > 1^177. It is found that the decay Z — > uuj gives a more 

f-*) . stringent bound on than that obtained from Z — > 1/1/77; the latter decay in 

turn gives limits on the neutrino-two-photon interaction that are less stringent than 
Q^' those obtained for a sterile neutrino I's from the analysis of f^iV i/sN conversion. 

Dh! 

The behavior recently observed of atmospheric [1] and solar [2] neutrinos pro- 
^ ! vides rather strong evidence that neutrinos have mass. This fact has renewed 

^ I the interest in neutrino electromagnetic properties, which have received con- 

siderable attention as they may shed light on some physics issues. In parti- 
cular, neutrino-one-photon interactions are of interest since they may play a 
key role in elucidating the solar neutrino puzzle, which can be explained in 
part by a large neutrino magnetic moment [3]. In the simplest extension of 
the standard model (SM), with the presence of massive neutrinos, one- loop 
radiative corrections induce a small magnetic moment proportional to the neu- 
trino mass m^, i.e. fx^ = 3 e G p ni,y / {8 \^ n'^) = 3 x 10~^^ rrii, hb [4], where 
my is to be expressed in KeV and /ib stands for the Bohr magneton. Several 
models have been advanced in order to induce neutrino magnetic moments 
as large as 10~^^-10~^° /xb [5], even with neutrino masses compatible with 
the mass square differences needed by atmospheric [1], solar [2] and the liquid 
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scintillation neutrino detector (LSND) [6] data. As for neutrino-two-photon in- 
teractions, they may have direct implications on several low- and high-energy 
reactions with astrophysical and cosmological interest [7] . For instance, a high 
annihilation rate of photons into a neutrino pair may explain the observed 
cooling of stars by neutrino emission [8] . In addition, there are other interest- 
ing processes involving ncutrino-two-photon interactions, such as — > 1/7, 
vv 77, and the neutrino double-radiative decay Vj ^'477. It is impor- 
tant to note that in the SM with massive neutrinos, the decay Vj — > 1/^77 is 
not severely suppressed by the GIM mechanism and can be the main decay 
channel as long as the Vj mass hes in the range of a few tenths of a MeV [9] . 

From the experimental side, the L3 collaboration searched for single-photon 
events near the Z pole at the CERN LEP collider and set a bound on the rare 
decay Z i/z/7 [10]. It was shown that the collected data impose a stringent 
constraint on the r neutrino magnetic moment [10-12]. In fact, the decay Z — > 
vT>^ can be a valuable tool to search for evidences of new physic since its rate is 
negligibly small in the SM [13]. By using the experimental bound on Z ^t^Ii 
an analysis in the framework of the effective Lagrangian approach (ELA) was 
carried out in Refs. [11,14] in order to constrain the operators that induce the 
couplings i^P7, 1/PZ7 and ZZ^. In regard to the rare decay Z — > vV^^, long 
ago the L3 and the OPAL collaborations looked for events with a lepton pair 
accompanied by a photon pair of large invariant mass [15]. After combining 
the data of both searches, the OPAL collaboration set an upper bound on the 
rate of Z — > vV^y. it was found that BR(Z 1/1/77) ^ 3.1 x 10~^. 

In the present letter we consider the possibility of obtaining indirect bounds 
on neutrino electromagnetic interactions from the experimental constraints on 
the decays Z — > uu'~^ and Z — > 1/1/77. Our main goal is to study these processes 
in a model independent way. We will also show that our results can be used 
to constrain the Vr magnetic moment and the decay Vj — > 1/^77. 

For the purpose of this letter we will consider the following effective interaction 

where /^^,. = /^i,.^,. is the Ui magnetic moment and [x^^^. {i ^ j) is the transi- 
tion magnetic moment. Although we will focus on Dirac neutrinos here, the 
discussion can be readily extended to Majorana neutrinos. 

As already mentioned, in Ref. [13] it was shown that the SM rate of the decay 

Z z/z/7 is unobservably small. Therefore, it represents an extraordinary 
mode to look for evidences of new physics arising from neutrino-one-photon 
interactions at a future e^e~ linear collider. The search for events with ener- 
getic single-photons along with missing energy at LEP was used by the L3 
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Fig. 1. Feynman diagrams contributing to the decay Z uu^ in the effective 
Lagrangian approach. The dots denote effective couplings. 

collaboration to set the bound BR(Z vD-f) < 10'^ [10]. Within the ELA, 

the rare decay Z — > z/z/7 can proceed through the Feynman diagrams shown 
in Fig. 1. For details of the analysis of these diagrams, we refer the reader 
to Refs. [11,14]. In particular, the experimental limit on Z — > 1/1/7 gives the 
following bound on the i>t magnetic moment 



This bound is in good agreement with that found by the L3 collaboration 
[10], and compares favorably with the bounds fi^,^ < 4 x 10~^ fi^ [16] and 
/ii,^ < 2.7 X 10~^ [17]. The former was obtained from low-energy experi- 
ments, whereas the latter was derived from the invisible width of the Z boson. 
Furthermore, our bound is close to the one obtained from a beam-dump ex- 
periment [18]. It is important to note that the most stringent bounds on the 
neutrino magnetic moment are obtained from chirality flip in supernova [19]. 

We now turn to examine the rare decay Z uP^^, which also can receive 
contributions from a neutrino-one-photon interaction through the Feynman 
diagrams depicted in Fig. 2. We would like to analyze if this decay is useful to 
bound the neutrino magnetic moment. Given the effective interaction of Eq. 
(1), the amplitude for the Feynman diagrams of Fig. 2 plus the crossed ones 
reads 



where the ellipsis stands for the crossed diagrams contribution, which can be 
obtained from the first term after the substitutions a ^ /5 and ki ^ k2. In 
principle, we must take into account all the neutrino species. However, in order 
to get an upper bound on /i,^^ , we will make a few assumptions for the sake of 
simplicity. First of all, we consider that the neutrino magnetic moment matrix 
is almost flavor diagonal, i.e. fi^,^ ^ fii^.i,. (j 7^ i). Secondly, we assume that the 
Ut- magnetic moment dominates over /ij,^ and /ij,^, i.e. there is the hierarchy 
3> Hu/j, 3> fJ^ue- III feet, the most stringent experimental bounds are fi^^ < 



1.1 X 10-^° hb [20], /i^^ < 7.4 X 10-9 i^B [20], and /x^, < 5.4 x 10"^ /xb [18]. 



Af^, < 2.62 X 10"^ AtB. 



(2) 



(3) 
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Bearing in mind these assumptions, the main contribution to the Z — > 
rate will arise from the Vr'^rl vertex. Therefore we can write M"^^^ as 

2 3 

M-^^ = ^ E Mp2)Tf'i^r{pi) h,k,„ (4) 

with 

rf^ = a-' + " a^^ - ^) " r Pl, (6) 
rf^ = a"^ + /^i) 7^ a'^^ (A + ^2) , (7) 

where we have neglected the z/^ mass; pi {P2), ^1,2 and p are the neutrino 
(antineutrino), photon and Z boson four-momenta; and Pl = {1 — 7^)/2 is 
the left-handed helicity projector. 

The transition amplitude can be squared by the usual trace technique. The 
result is too lengthy to be shown here. The squared amplitude can then be 
integrated over the four-body phase space with the aid of the Monte Carlo in- 
tegration method [21]. In order to cross-check our results, we used two different 
methods for the evaluation of the Z — > 1/^77 decay rate. In the first method we 
squared the amplitude and then used a Monte Carlo event generator to carry 
out the numerical integration [22]. As far as the second method is concerned, 
we implemented the 1/1/7 and z/i/77 interactions into the CALCHEP program 
[23], which automatically generates the respective set of Feynman diagrams, 
squares the matrix elements, and integrates over the phase space. There was 
nice agreement between the results obtained by both of these methods. 

Under the assumptions discussed above, we can obtain the following estimate 
for the Z — > vPj^ rate 

BR(Z ^ i/P77) = 1.749 x 10^^ ( , (8) 

which along with the experimental bound BR(Z vV^^) < 3.1 x 10~^ yield 

li,^ < 6.488 X 10-^ /XB. (9) 

which is just one order of magnitud below than the bound obtained from the 
three body decay Z — > vPj [cf. Eq. (2)]. The rare decay Z — > 1/P77 is however 



4 




Fig. 2. Feynman diagrams contributing to the decay Z 1/1/77. The crossed dia- 
grams are not shown. The dots denote effective couphngs. 

more sensitive to the value of the neutrino magnetic moment. In this respect, 
it is interesting if we take a different approach and use the most stringent 
experimental bound on yU^^ [18] to constrain the rare decays Z — > 1/P7 and 
Z i/i/jj, in which case we are led to 

BR{Z uP-f) < 6.917 X 10~^ (10) 

BR{Z I/P77) < 1.487 x 10"^^ (11) 

Before proceeding, we would hke to stress that the procedure described above 
can be employed to get bounds on the r neutrino transition magnetic moments 
f^urUi {i = 6) A^)) which also have been the source of interest [24]. In that case, 
we would have to include all the ^vrVi contributions into the Z vv^ and 
Z z/z/77 rates, and the upper bound on each \ivrVi would be obtained after 
dropping the remaining contributions. Since the results are insensitive to the 
neutrino mass, provided that rrii, <S m^, it follows that the bounds of Eqs. (2) 
and (9) also apply to [1^^^^. Of course the same is true for any /Xj^.j^ , but the 
bounds are very weak for and v^^ as compared to other results appearing 
in the literature. 

Now we will analyze the impact of the z/z/77 coupling on the rare decay 
Z — > z/z/77. While the neutrino-one-photon interaction is strictly vanishing 
for massless neutrinos, the ncutrino-two-photon interaction can have a nonzero 
value even if neutrinos are massless. Long ago, it was shown that if massless 
neutrinos interact locally with charged leptons the neutrino- two-photon vertex 
vanishes [25]. It is possible that neutrinos are kept massless but interact di- 
rectly with gauge bosons, as in the SM. In that case there is indeed a nonzero 
neutrino-two-photon coupling, which arises at the one-loop level and is neg- 
ligibly small, of order 0{Gy) [26]. Nevertheless, as pointed out in Ref. [27], 
the introduction of massive neutrinos can enhance dramatically the vv^ and 
i/p77 vertices. While in several SM extensions the neutrino-two-photon inter- 
action is proportional to the neutrino mass, there are some extensions, such 
as left-right symmetric theories and the Zee model, where it is proportional 
to a heavy Higgs scalar mass [28]. This fact may give rise to a significant en- 
hancement of the i/p77 coupling, which can be parametrized in the following 



5 



way at the lowest dimension [9,27,28] 

^^^-^T7 = 4^ {^lPl + cxIPr) ujF^T^,, (12) 

where are dimensionless couphng constants and A is the new physics 
scale. Given this interaction, the procedure described before can be used to 
obtain the contribution to the decay Z — ^ z/z/77 from diagram 2(d) plus that in 
which the photon pair emerges from the neutrino. We can write the respective 
transition amplitude as follows 

^''^'^ = T^3''i(P2)^'^AP^) ^^""^ ^lA^2p, (13) 

with 

= at Pl - A) Y + Oil Pr r - ^2) ■ (14) 

Again, we have neglected the neutrino masses since the result is insensitive to 
them. The squared amplitude can be written in a very short way: 

\M? = ly^t'ji {K?nP.) + WM'FM) . (15) 

with 

F(P.) = 7^ U -Hp -P.?- '-^^^ (P - ■ P.) . (16) 

(P - Pi) V Pi • P2 J 

From the last expressions, the Z — z/z/77 decay rate can be obtained after 
Monte Carlo integration. Hereafter, we will consider the contributions from 
the three SM neutrino species. The resulting branching fraction is thus given 

by 

BR{Z ^ i/P77) = 1.092 x 10^ EE (l^Lf + I^^T) 

i j 

where the sums run over z/g, z/^, and Vt- After using the experimental bound 
on Z — > i/p77, we are left with 

'EE(l«^r+l«^r)^2.85xlO-^ (18) 

i j 
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Fig. 3. Invariant mass distribution of the photon pair {X^^ = \J + h^)^ jmz) in 
the rare decay Z — > vv^^. The contributions from the vertices vv^ and 1/1/77 ^"^^ 
shown separately. 

This bound is weaker than that obtained for a sterile neutrino Vg from the 
analysis of the PrimakofT effect on the process of v^N — > conversion in 
the external Coulomb field of the nucleus N [29]. 

At this point, we would like to note some interesting features of the photon 
energy and invariant mass distributions of the decay Z vv^^. In Fig 3 
we have plotted the distribution of the invariant mass of the photon pair 
when the contribution from either vertex vv^ or vv^^ is considered at a time. 
When only the 1/P7 interaction contributes, the invariant mass peaks around 
X^-y = 1/4; on the other hand, when the vv^^ vertex alone contributes, the 
peak is located around = 1/2. A similar situation is observed in Fig. 
3, where we have plotted the energy distribution of the photon pair, and 
in Fig. 5, where it is shown the energy distribution of an isolated photon. 
We can observe that in both plots the peak of the curve accounting for the 
vv^ contribution is shifted to the left with respect to the curve resulting 
from the z/z/77 contribution. Therefore, in principle a proper set of cuts would 
allow us to distinguish between the contributions from each vertex. A more 
comprehensive analysis is however beyond the present letter. 

Finally, we will show that the decay Z — > 1/1/77 can also be used to bound the 
neutrino double radiative decay Vj ^'^77. Neglecting the z/j mass and after 
some calculation one can obtain the following expression for the Vj — > vf^^ 
decay width 
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Fig. 4. Energy distribution (Ky^ = (Ei + E2)/m,z) of the photon pair in the rare 
decay Z — > The contributions from the vertices uv'y and 1^1/77 are shown 

separately. 
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Fig. 5. Energy distribution {Y^ = 2Ei/mz) of a single photon in the rare decay 
Z — > vu^^. The contributions from the vertices vi>^ and 1^2/77 are shown separately. 



which yields 



= 1.59 X 10 



-3 



(I 



a 



ij\2 



+ \a 



ij|2i 

R\ ; 



"iGeV 
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6 r 



-1 



IMeV 



(20) 



The constraint of Eq. (18) can then be translated into a lower bound on the 
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Uj lifetime 

T^. > 1.79 X 10 

which is true provided that rrii,. lies in the range of a few tenths of a MeV, 
since in that case Uj 1/^77 is the dominant decay channel [9] . Our bound is 
one order of magnitude below than the one previously found for the lifetime 
of a sterile neutrino [29] . 

In closing we emphasize that the rare decay Z 1/P7 gives rise to a bound on 
the magnetic moment that is in excellent agreement with other ones found 
recently. The bound obtained from the rare decay Z uu'-f-f is just one order 
of magnitude below. It must be stressed that the current experimental bound 
on the Z — >■ 1/1/77 branching ratio is somewhat weak. In fact, it is of the same 
order of magnitude than that on Z — > vP^. Some improvement is expected 
from the data to be collected at a future linear collider. One important feature 
of the decay Z uu'-f'-f is that it can also be used to bound the ncutrino- 
two-photon interaction. In this respect, the resulting bound is weaker than 
that derived for a sterile neutrino from the process of ly^N — > UgN conversion 
in the external Coulomb field of the nucleus N. Finally, the bound on the 
neutrino-two-photon interaction also allowed us to constrain the width of the 
decay Uj i^i'yy, which in turn can be translated into a constrain on the uj 
lifetime as long as m„. lies in the range of a few hundreds of KeV. The main 
advantage of our procedure is that it is model-independent and relies on a few 
assumptions. 
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